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Previous studies have demonstrated tliat tlie difference between resting-state brain 
activations depends on winether the subject was eyes open (EO) or eyes closed 
(EC). However, whether the spontaneous fluctuations are directly related to these 
two different resting states are still largely unclear. In the present study, we acquired 
resting-state functional magnetic resonance imaging data from 24 healthy subjects 
(11 males, 20.17 ± 2.74 years) under the EO and EC states. The amplitude of the 
spontaneous brain activity in low-frequency band was subsequently investigated by using 
the metric of fractional amplitude of low frequency fluctuation (fALFF) for each subject 
under each state. A support vector machine (SVM) analysis was then applied to evaluate 
whether the category of resting states could be determined from the brain spontaneous 
fluctuations. We demonstrated that these two resting states could be decoded from 
the identified pattern of brain spontaneous fluctuations, predominantly based on fALFF 
in the sensorimotor module. Specifically, we observed prominent relationships between 
increased fALFF for EC and decreased fALFF for EO in sensorimotor regions. Overall, the 
present results indicate that a SVM performs well in the discrimination between the brain 
spontaneous fluctuations of distinct resting states and provide new insight into the neural 
substrate of the resting states during EC and EO. 

Keywords: resting-state fMRI, fractional amplitude of low-frequency fluctuation (fALFF), support vector machine 
(SVM), eyes closed, eyes open 



INTRODUCTION 

The modulation of resting-state brain activity has attracted great 
attention in recent years (Htifner et al., 2009), and previous 
studies have demonstrated that resting-state brain activity is dif- 
ferentially affected based on whether the subjects have their eyes 
open (EO) or eyes closed (EC). Until recently, the widespread 
influences of EO versus EC on resting-state activity have been 
demonstrated using methods ranging from regional spontaneous 
activity (Marx et al., 2003, 2004; Mcavoy et al, 2008; Bianciardi 
et al., 2009) to functional connectivity (Yan et al., 2009; Zou et al., 
2009; Wu et al., 2010; Mcavoy et al, 2012) and to network topo- 
logical organization (Jao et al., 2013; Xu et al., 2014). Despite these 
advances, little is known about whether there are spontaneous 
fluctuations that are directly related to eye behavior states. 

Several functional magnetic resonance imaging (fMRI) studies 
have already shown that spontaneous brain oscillations are dif- 
ferent between resting states with EC and EO in many regions, 
including the visual cortex (Raichle et al., 2001; Uludag et al., 
2004), auditory cortex (Marx et al., 2003; Qin et al., 2013), 
somatosensory cortex, fronto-parietal attentional regions, and 



default mode network (DMN) (Yan et al, 2009). Importantly, 
alterations of spontaneous activity associated with different eye 
behavior states vary between EC and EO. For example, our recent 
study (Xu et al., 2014) found that the attentional, ocular motor, 
and arousal systems showed higher regional nodal properties 
(e.g., efficiency, degree, and betweenness centrality) in EO than in 
EC; in contrast, the visual, auditory, and somatosensory systems 
and parts of the DMN (Raichle et al., 2001) showed significantly 
lower regional nodal properties in EO than in EC. These find- 
ings may confirm the existence of an "exteroceptive" mental state, 
characterized by attentional and ocular motor activity during EO, 
and an "interoceptive" mental state, characterized by imagination 
and multisensory activity during EC (Marx et al., 2004). It should 
be noted that the notion of exteroceptive and interoceptive men- 
tal states is derived from the direct measurement of spontaneous 
brain oscillations (Marx et al., 2004). In addition, alterations in 
spontaneous activity (e.g., amplitude of low frequency fluctu- 
ation, ALFF; fractional ALFF, fALFF) between the EO and EC 
states have been found in the visual cortex, paracentral lobule 
(PCL) (Yang et al, 2007), middle temporal gyrus, anterior insula. 



Frontiers in Human Neuroscience 



www.frontlersin.org 



August 2014 | Volume 8 | Article 645 | 1 



Liang et al. 



Sensorimotor fALFF of eyes conditions 



and DMN regions (Yan et al, 2009). Although the widely dis- 
tributed differences of brain activity are related to these two 
mental states, the knowledge whether the spontaneous fluctua- 
tions in these brain regions could combine to be directly related 
to the two states is still largely limited. Thereby, a sophisticated 
approach is consequently necessary to assess the link of the spon- 
taneous fluctuations to the eye behavior states, which will pro- 
mote the investigation of fundamental neural differences between 
EC and EO. 

Recently, machine learning approach has been widely applied 
to investigate the relationship between the distributed functional 
signal and the brain states (Craddock et al., 2009; Poldrack et al., 
2009; Dosenbach et al., 2010; Shen et al, 2010; Naselaris et al, 
2011; Zeng et al, 2012; Vergun et al, 2013). By building a math- 
ematical model (e.g., support vector machine, SVM), we could 
predict the brain states from the spatially distributed sponta- 
neous activities and further assess their relationship directly. In 
the present study, we sought to identify whether the spontaneous 
fluctuation amplitude directly link to the EC and EO by using the 
SVM approach. Toward this end, we recruited 24 healthy subjects 
and collected their resting state fMRI data. By using fALFF, we 
first localized those brain regions showing differences in sponta- 
neous resting-state activity between EC and EO. Next, a machine 
learning method was further applied to examine whether there 
exists a pattern of the spontaneous fluctuation amplitude which 
is directly related to the eye behavior states. 

MATERIALS AND METHODS 
SUBJECTS 

A total of 24 right-handed healthy undergraduates/postgraduates 
(11 males/13 females, 20.17 ± 2.74 years) were recruited for 
the present study. No subjects had a history of neurological or 
psychiatric disorders or head injury. This study was approved 
by the Institutional Review Board of Beijing Normal University 
Imaging Center for Brain Research. Written informed consent 
was obtained from each subject before the experiment. 

DATA ACQUISITION 

All MRI data were acquired on a 3T Siemens Trio TIM MR 
scanner powered with a total imaging matrix (TIM) tech- 
nique at the Imaging Center for Brain Research, Beijing Normal 
University. The resting-state functional magnetic resonance imag- 
ing (rsfMRI) data were acquired using a 12-channel phased array 
receiver-only head coil. The rsfMRI data were obtained using 
a gradient-echo EPI sequence with the following two param- 
eter sets: (A) TR = 2000 ms, TE = 30 ms, 33 transverse slices, 
slice thickness = 3.5 mm, gap = 0.7 mm, flip angle = 90°, 
FOV = 224 mm x 224 mm, matrix = 64 x 64, and 240 volumes 
covering the whole brain; (B) TR = 3000 ms, 40 transverse slices, 
slice thickness = 3.5 mm, no gap, 160 volumes, and the other 
parameters identical to those of (A). For each parameter set, all 
subjects underwent the rsfMRI scan under both conditions (EC 
and EO) in turn. Thus, there were four types of rsfMRI scan: 
AO, AC, BO, and BC. For each subject, all four scans were per- 
formed in the same session. The order of data acquisition was 
counterbalanced across all subjects. In addition, we also acquired 
high-resolution 3D brain structural images for each subject by 



using MP-RAGE sequence with the implementation of parallel 
imaging scheme GRAPPA (GeneRalized Autocalibrating Partially 
Parallel Acquisitions) (Griswold et al., 2002) and the acceleration 
factor 2. 

DATA PREPROCESSING 

For the present study, we analyzed the rsfMRI data correspond- 
ing to TR = 2000 ms only. Data preprocessing was performed 
using DPARSF (http://www.restfmri.net/forum/DPARSF) with 
SPM8 (http://www.fil.ion.ucl.ac.uk/spm/). We discarded the first 
10 volumes of the functional images to account for signal equilib- 
rium and the subjects' adaptation to the scanning environment. 
We then corrected the remaining functional images for within- 
scan acquisition time differences between slices and realigned 
all images to the first volume to correct for head motion. This 
realignment procedure provided us with a record of head motion 
during each rsfMRI scan. None of the subjects were excluded for 
excessive head motion based on criteria of >2mm displacement 
or an angular rotation of >2° in any direction. The summary 
scalars of both gross (maximum and root mean square) and micro 
(mean frame-wise displacement) head motion were matched 
between the two conditions (all p > 0.14). The corrected func- 
tional images were subsequently spatially normalized to the MNI 
standard template using an optimum 12-parameter affine trans- 
formation and non-linear deformations and were re-sampled to a 
voxel size of 3 x 3 x 3 mm^ . Finally, the linear trend of the fMRI 
data was estimated with a least-square fitting of a straight line and 
was removed from the functional data. 

CALCULATION OF fALFF 

We calculated the fALFF in a voxel-wise manner following the 
procedures described in previous studies (Zang et al., 2007; Zou 
et al., 2008; Song et al., 201 1). After data preprocessing, we trans- 
formed the data time series to the frequency domain using fast 
Fourier transformation (EFT), obtained the power spectrum, and 
then averaged across 0.01-0.08 Hz at each voxel. The power of a 
given frequency is proportional to the square of the amplitude of 
this frequency component in the original time series in the time 
domain. This averaged square root was defined as the ALFF. Next, 
the fALFF was computed as the ratio of the ALFF in the range 
of 0.01-0.08 Hz to that of the entire frequency range (0-0.25 Hz) 
for each voxel. The fALFF map was smoothed (fuU width at half 
maximum = 6 mm) and the normalized fALFF map was obtained 
by dividing fALFF at each voxel by each individual's global mean 
fALFF value. 

INDIVIDUAL fALFF MAPS 

Using an existing template that was functionally defined based 
on a meta-analysis in a previous study (Dosenbach et al., 2010), 
we defined 160 cortical regions of interest (ROIs) for the entire 
cortex (Table SI). This template broadly covers most of the cor- 
tex and cerebellum and has been used in several studies, such 
as brain functional network studies (Wang et al, 2011; Hwang 
et al, 2013). We first extracted the fALFF value at each voxel for 
each resting condition (EO or EC). We then calculated the mean 
fALFF for each ROI by averaging the fALFF values across aU vox- 
els for the selected ROI in each subject. Finally, we divided the 160 
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ROIs into six modules (default network, fronto-parietal, cingulo- 
opercular, sensorimotor, occipital, and cerebellum) according to a 
meta-analysis of cognitive experimental fMRI studies (Dosenbach 
et al., 2010). Thus, individual fALFF maps were constructed, and 
fALFF values for the 160 ROIs were obtained for each resting 
condition. 

MACHINE LEARNING METHOD 

The LIBSVM toolbox (http://www.csie.ntu.edu.tw/~cjlin/ 
libsvm/) with the linear support vector classification (SVC) was 
applied as the classifier for the multivariate pattern analysis 
(MVPA) (Zhang et al., 2013). The flowchart of the analysis 
stream for the MVPA method is shown in Figure 1 . 

The first step was feature selection. To reduce the dimension- 
ality of the discriminatory patterns, we selected features using 
a nonparametric two-sample Wilcoxon signed-rank test (Chen 
et al., 2011). First, we calculated the mean fALFF for each ROI 
by averaging across all subjects in two resting conditions (EO and 
EC) separately. Then, the differences among all the features, rep- 
resented by the mean fALFF value for each ROI (a total of 160 
features) between the two resting conditions (EO and EC), were 
estimated to have a normal distribution and were selected as the 
discriminatory features for MVPA. 

To extract highly discriminatory brain areas, we used four 
rankings of the difference scores (Z-scores) of the Wilcoxon 
signed-rank test as the criteria for the feature selection: |Z| > 1.96 



(-0.05 < p < 0.05), Z > 1.96 (p < 0.05), Z > 2.25 (p < 0.01), 
and Z < —1.96 (p > —0.05). Using these four criteria, we were 
able to systematically decrease the number of different regions in 
the fALFF maps between the two conditions. The selected features 
were the same for all subjects based on different Z criteria. 

The second step was pattern assembly. Four pattern types were 
identified from the selected discriminatory features according to 
different criteria. For each discriminatory pattern, we ranked the 
features of each subject for each resting condition (EO or EC) in 
a row vector to constitute the sample data, and individual rows 
were combined into a sample dataset. For each of the 24 subjects, 
the rsfMRI scan was performed under two resting conditions, EC 
and EO, resulting in a dataset containing 48 samples. Finally, the 
samples for each resting condition were randomly divided into 
two sets of 12 samples each. Half of the data was used to construct 
the training dataset (24 samples across two conditions), and the 
other half was used as the testing dataset (24 samples across two 
conditions). The aim was to improve the reliability of the classifier 
validation as much as possible. 

The third step was classifier verification. The sample data 
in the training and testing datasets were labeled with category 
labels (EC = 1 and EO = — 1). The SVC classifier was trained 
on the training dataset and the corresponding category labels, 
and the performance of the trained classifier was tested using 
the testing dataset and the corresponding category labels. Using 
a receiver operating characteristic (ROC) curve, we evaluated the 




FIGURE 1 I Flowchart of MVPA. The fALFF values of 160 ROIs were 
used as the features. A nonparametric two-sample Wilcoxon signed-rank 
test was performed to detect discriminatory patterns between two 
conditions (EC and EO). The selected features from each of the 12 
subjects for each condition were ranked as a sample dataset, which 
resulted in two sample datasets for each condition: one was used to 



construct the training dataset (24 data samples across two conditions) 
and another one for the testing dataset (24 data samples across two 
conditions). The SVM classifier was trained with the training dataset, and 
the performance of the trained classifier was tested with the testing 
dataset. In step 3, stars and triangles represent the samples of the two 
conditions. 
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performance of the trained classifier based on its accuracy which 
was estimated as the percentage of the correct predicted category 
label. 

RESULTS 

fALFF IN EO AND EC 

There was a significant difference in the mean fALFF between the 
two resting conditions, EC and EO, using a two-sample Wilcoxon 
signed-rank test (|Z| > 1.96; Figure2). Compared with EO, the 
fALFF values in EC were significantly increased in the senso- 
rimotor module and cingulo-opercular region but significantly 
decreased in the fronto-parietal cortex, occipital cortex, and cere- 
bellum. Additionally, we noted that fALFF changes between EO 
and EC were far from consistent in the DMN regions. Compared 
with EO, the fALFF value in EC was higher in the left inferior pari- 
etal supramarginal gyrus (IPS. L), but lower in the left occipital 
cortex and ventromedial prefrontal cortex (vmPFC.L) (Table 1). 

SVM CLASSIFIER PERFORMANCE 

We found that the eye conditions could be decoded from the brain 
fALFF map using a SVM composed of 28 discriminatory brain 
regions(|Z| > 1.96) (Figures). Theaccuracyofthisclassifierwas 
85% (Table 1). 

Using a subset of brain regions, the accuracy of the classifier 
was 84% (21 brain regions with Z > 1.96) and 88% (13 brain 
regions with Z > 2.25). However, the classifier showed poor per- 
formance (accuracy = 63%) based on a pattern of 7 brain regions 
(Z < —1.96), i.e., the fALFF of EC was significantly lower than 
that of EO. 

Figure 4 shows the ROC curve of the neural activity patterns 
associated with different brain regions (|Z| > 1.96, Z > 1.96, 
Z < —1.96, and Z > 2.25). This figure clearly demonstrates that 
the pattern of brain regions with Z > 1.96 was similar to the 
pattern of regions with |Z| > 1 .96 in decoding the eye conditions. 

For Z > 1.96, the brain regions were widely distributed across 
the sensorimotor module, DMN, and cingulo-opercular region, 
with nearly 86% of these brain regions located in the senso- 
rimotor module. For Z < —1.96, however, the brain regions 
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FIGURE 2 1 Brain regions involved in the discriminative pattern 


corresponding to |Z| > 1.96. The details of the brain regions and modules 


are listed in Table SI 









were mainly located in the DMN, fronto-parietal cortex, occipital 
cortex, and cerebellum. 

DISCUSSION 

The present study investigated whether the fALFF could be 
directly related to distinct resting conditions (i.e., EC and EO) 
using a SVM. The main findings can be summarized as fol- 
lows: (i) compared with EO, the EC condition induced higher 
fALFF values primarily in the sensorimotor cortex and cingulo- 
opercular region but lower fALFF values in the fronto-parietal 
cortex, occipital cortex, and cerebellum; (ii) the eye conditions 
could be decoded from the fALFF of the sensorimotor regions 
using a SVM; and (iii) increased or decreased fALFF in the senso- 
rimotor cortex was directly related to the EC and EO conditions, 
respectively. 

fALFF OF DIFFERENT RESTING STATES 

The present study demonstrated that changes in fALFF were 
spatially distributed and predominantly located in the poste- 
rior occipital cortex, DMN-associated regions (e.g., intraparietal 
sulcus [IPS]), fronto-parietal regions, sensorimotor cortex, cere- 
bellum, basal ganglia, and middle insula (Table 1). These results 
are consistent with previous studies (Yang et al, 2007; Yan et al., 
2009; Jao et al, 2013), which indicated that the ALFF can reflect 
different cerebral physiological states (Yang et al., 2007; Zang 
et al., 2007). Furthermore, differences in the ALFF between the 
EC and EO states have previously been identified in regions such 
as the right PCL, visual cortex (Yang et al., 2007), and regions 
associated with the DMN (e.g., posterior cingulate cortex [PCC] 
and medial prefrontal cortex [MPFC]; (Yan et al, 2009). 

Deactivation of the DMN is well known to account for a sub- 
stantial portion of resting-state brain activity (Raichle et al., 2001; 
Fox et al., 2005; Fransson, 2005), and researchers have shown that 
neural activities in the DMN, one aspect of functional connectiv- 
ity, were similar across different resting states such as EC and EO 
(Fox et al, 2005; Fransson, 2005). Accordingly, Yang et al. (2007) 
found no significant differences in the ALFF between the EC and 
EO states in the PCC, which is a key region of the DMN. These 
findings supported the conclusion that the DMN was not easily 
disrupted by visual information compared with the visual cortex 
(Greicius et al, 2003). 

However, there have been discrepant findings regarding the 
neural activity of the DMN. For example, Yan et al. (2009) found 
that EO induced higher functional connectivity in the DMN and 
increased ALFF compared with EC. This observation suggests that 
increases in spontaneous neural activity may be related to greater 
involvement of non-specific neural activity during EO (e.g., non- 
goal-directed visual information gathering and evaluation or 
mind-wandering and daydreaming). Using nonparametric meth- 
ods, our present data indicate that changes in the activity of 
DMN regions were inconsistent between the EC and EO states. 
We found that the IPS showed higher fALFF during EC whereas 
some parts of the occipital cortex showed lower fALFF during 
EC (Table 1). This finding may suggest the existence of differen- 
tial responses of different DMN regions between the EC and EO 
states, and the identification of these DMN changes may depend 
on specific statistical methods. 
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Table 1 | Brain regions involved in the discriminatory pattern with |Z| > 1.96. 
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The (x, y, z) coordinates correspond to the peak voxel in MNI space. The Z-vaiue indicates the difference scores, which were calculated from the nonparametric 
two-sample Wilcoxon signed-rank test, and a positive value indicates that the Z-value of EC was higher than that of EO (and vice versa). Six modules were evaluated, 
including default, fronto-pahetal, cingulo-opercular, sensorimotor, occipital, and cerebellum. L (R), left (hght) hemisphere. 



Similarly, there are also contradictory views regarding 
spontaneous neural activity in the visual cortex between the EC 
and EO states. For example, a lower ALFF in the bilateral visual 
cortex and higher ALFF in the right PCL were observed in EC 
compared with EO (Yang et al., 2007). One possible interpreta- 
tion is that the visual cortex was influenced by visual input during 
EO (Raichle et al, 2001; Uludag et al, 2004). In contrast, several 
studies indicated that EC induced greater changes in the sponta- 
neous neural activity in visual regions (Marx et al., 2003, 2004). In 
addition to these local changes in brain activity, lao et al. (2013) 
found that the mean fALFF of the whole brain was significantly 
higher during EC than during EO. 

Considering the widely distributed differences in fALFF 
between EC and EO and conflicting findings of neural activity 
changes within the same region, we speculate that the spe- 
cific pattern of results may depend on the statistical methods. 
Although several statistical methods (e.g., univariate, bivariate, 
and multivariate methods) have already been used to detect 
changes in the fALFF between eye conditions (Yang et al, 2007; 



Yan et al, 2009; Mcavoy et al, 2012; Jao et al, 2013), to the best 
of our knowledge, few studies have used multivariate methods to 
identify and decode discriminatory fALFF patterns with regard 
to different eye conditions. Of note, considering not all of the 
changed brain activity regions were directly related to the eyes 
behavior states, it is necessary to investigate which characteris- 
tic changes of the fALFF allow for the accurate discrimination 
between the EC and EO states. 

fALFF PAHERNS OF EYE CONDITIONS 

The advantages of multivariate methods have been demon- 
strated in many previous studies (for review, see Norman et al, 
2006). Using SVM methods, neural activity patterns have been 
identified in fMRI data that allow for the classification of dif- 
ferent cognitive tasks (Poldrack et al., 2009), various diseases 
(Craddock et al, 2009; Poldrack et al, 2009; Shen et al., 2010), 
and brain aging (Vergun et al, 2013). Inspired by these previ- 
ous observations, we used a multivariate method, SVM, based 
on a nonparametric test approach (Chen et al., 2011), to detect 
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FIGURE 3 I ROC curves of the performance of SVM classifiers corresponding to different neural fALFF patterns. Tine yellow line corresponds to the 
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FIGURE 4 I Mean fALFF values of each brain region in the neural fALFF 
map with |Z| > 1.96 related to the EC and EO states. fALFF was 
calculated as the ratio of the ALFF value in the range of 0.01-0.08 Hz to that 
of the entire frequency range (0-0.25 Hz). The brain regions and 
corresponding modules are same to those in Figure 2. 



the relationship between the whole brain fALFF and the eye con- 
ditions (for details, see "Materials and methods"). The present 
study found that differences in the fALFF between the EC and 
EO states were observed in regions associated with the DMN 
but that the pattern of changes across these two conditions was 
inconsistent for this region (Table 1). In addition to the visual 
cortex and DMN regions, we also observed widespread regions. 



including the sensorimotor cortex, basal ganglia, middle insula, 
fronto-parietal cortex, and cerebellum, with distinct fALFF val- 
ues between the EO and EC states (Figure 3). Our results are 
consistent with those of previous studies, in which a number of 
distributed regions were affected by eye conditions. For exam- 
ple, EO activated ocular motor and attentional systems associated 
with "exteroceptive" mental states (Marx et al, 2004), whereas 
EC activated the visual, somatosensory, and auditory systems 
and caused significantly higher fALFF in posterior brain regions 
associated with "interoceptive" mental states (Marx et al., 2004). 

Although changes in the fALFF were observed in many regions, 
not all the changes could be used to discriminate between the 
eye conditions. In the present study, we found that the fALFF 
in the sensorimotor cortex was a significant contributor to the 
whole-brain fALFF map of EC compared with EO. This finding is 
consistent with a previous study, which showed that brain regions 
associated with sensory processing were easily influenced by the 
eye behavior states (Marx et al., 2003). 

Moreover, we found that the classifier performance was similar 
for the fALFF patterns corresponding to Z > 1.96 and |Z| > 1.96 
but that the classifier performance was much lower for the fALFF 
pattern corresponding to Z < —1.96. We found that sensorimo- 
tor regions composed approximately 86% of all regions in the 
fALFF pattern of Z > 1.96. The statistical analysis indicated that 
fALFF was significantly higher in sensorimotor regions in the 
EC state than in the EO state (Figure 4). There exists a possi- 
bility that EO leads to a suppression of sensory modalities to 
allocate resources to exteroceptive processes (Xu et al., 2014). 
In the current study, we further suggest that this neurological 
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process specific to sensorimotor regions may play a vital role in 
the predicting of mental states. 

EYE CONDITIONS AND MENTAL STATES 

Although exteroceptive and interoceptive states might not be 
directly manipulated by eye conditions, the current results and 
previous findings have demonstrated that EO and EC correspond 
to these two states and that each is associated with distinct brain 
activity patterns (Marx et al., 2003, 2004). Consistent with stud- 
ies that involved interoceptive manipulations (Critchley et al., 
2004; Simmons et al, 2013), sensorimotor regions including 
the insula, fi-onto-parietal regions, and precentral gyrus showed 
higher spontaneous activity during EC than during EO (Table 1). 
Activity in sensorimotor regions was previously observed to be 
associated with mental imagery during a resting state (Mazard 
et al, 2002, 2005; Assaf et al, 2006). Wang et al. (2008) found 
that spontaneous resting activity of the sensorimotor cortex was 
associated with that of the primary visual cortex when subjects 
remained in the EC state, suggesting that mental imagery occurs 
in the resting brain in the absence of tasks. Because sensorimo- 
tor activity patterns can accurately differentiate between these 
two states, the results of the present study suggest that there may 
be a unique mental state associated with the resting state dur- 
ing EC compared with task-dependent conditions and the EO 
resting state. The neural ensemble of the resting state during EC 
may provide an approach to explore the neural mechanisms of 
interoceptive states. 

LIMITATIONS 

This study contains several limitations that should be addressed 
in future studies. First, we concentrated on the amplitude of 
brain spontaneous fluctuations, and explored its association with 
two resting state conditions (i.e., EC and EO) with the help 
of the SVM in the present study. Other dimensions of spon- 
taneous activity, such as functional connectivity and network 
topological organization of spontaneous activity, should also be 
considered in future studies. Second, the whole-brain fALFF map 
was constructed based on 160 ROIs derived from a previous 
meta-analysis. It may be an interesting topic in future research 
to investigate whether the findings of the present study were 
dependent on the spatial resolution of the fALFF (e.g., voxel- 
wise, and AAL-90 or AAL-1024). Third, for the reason that 
the EC and EO conditions provide EEC measures differing in 
topography as well as power levels (Barry et al., 2007, 2009), 
the influence of the frequency band should be further consid- 
ered in the investigation of the neural substrate underlying EC 
and EO. 

CONCLUSION 

In summary, the present study investigated the associations 
between the amplitude of brain spontaneous fluctuations and the 
EO and EC conditions using a machine learning method. We 
found a direct relationship between different eye conditions and 
the fALFF during a resting state, in which most of the regions were 
located in the sensorimotor cortex. The fALFF changes (increases 
and decreases) in these sensorimotor regions were significantly 
modulated by the eye behavior states. Our findings may provide 



new insights into the neural substrate of distinct resting state 
activity associated with different eye states. 

ACKNOWLEDGMENTS 

This work was supported by the funding of National Natural 
Science Foundation of China under Grant numbers: 31371049, 
81071149, 81030028, and 81271548; the National Social Science 
Foundation of China under grant number BBAl 10411; the 
Provincial Department of Construction Project of Key Research 
Base of Humanities and Social Sciences of China under grant 
number 11JJD190003; and the Scientific Research Foundation of 
Graduate School of South China Normal University. 

REFERENCES 

Assaf, M., Calhoun, V. D., Kuzu, C. H., Kraut, M. A., Rivkin, P. R., Hart, J. Jr., 

et al. (2006). Neural correlates of the object-recall process in semantic memory. 

Psychiatry Res. 147, 115-126. doi: 10.1016/j.pscychresns.2006.01.002 
Barry, R. J., Clarke, A. R., Johnstone, S. J., and Brown, C. R. (2009). EEC differ- 
ences in children between eyes-closed and eyes-open resting conditions. Clin. 

Neurophysiol. 120, 1806-1811. doi: 10.1016/j.clinph.2009.08.006 
Barry R. J., Clarke, A. R., Johnstone, S. J., Magee, C. A., and Rushby, J. A. (2007). 

EEC differences between eyes-closed and eyes-open resting conditions. Clin. 

Neurophysiol 118,2765-2773. doi: 10.1016/j.clinph.2007.07.028 
Bianciardi, M., Fukunaga, JVI., Van Gelderen, R, Horovitz, S. G., De Zwart, 

J. A., and Duyn, J. H. (2009). Modulation of spontaneous fMRI activity 

in human visual cortex by behavioral state. Neuroimage 45, 160-168. doi: 

10.1016/j.neuroimage.2008.10.034 
Chen, G., Ward, B. D., Xie, C, Li, W., Wu, Z., Jones, J. L., et al. (201 1). Classification 

of Alzheimer disease, mild cognitive impairment, and normal cognitive status 

with large-scale network analysis based on resting-state functional MR imaging. 

Radiology 259, 213-221. doi: 10.1148/radioL10100734 
Craddock, R. C, Holtzheimer, R E. 3rd., Hu, X. R, and Mayberg, H. S. (2009). 

Disease state prediction from resting state functional connectivity. Magn. Reson. 

Med. 62, 1619-1628. doi: 10.1002/mrm.22159 
Critchley, H. D., Wiens, S., Rotshtein, R, Ohman, A., and Dolan, R.J. (2004). 

Neural systems supporting interoceptive awareness. Nat. Neurosci. 7, 189-195. 

doi: 10.1038/nnll76 

Dosenbach, N. U., Nardos, B., Cohen, A. L., Fair, D. A., Power, J. D., Church, J. A., 
et al. (2010). Prediction of individual brain maturity using fMRI. Science 329, 
1358-1361. doi: 10.1126/science.ll94144 

Fox, M. D., Snyder, A. Z., Vincent, J. L., Corbetta, M., Van Essen, D. C, and Raichle, 
M. E. (2005). The human brain is intrinsically organized into dynamic, anticor- 
related functional networks. Proc. Natl. Acad. Sci. U.S.A. 102, 9673-9678. doi: 
10.1073/pnas.0504136102 

Fransson, P. (2005). Spontaneous low-frequency BOLD signal fluctuations: an 
fMRI investigation of the resting-state default mode of brain function hypothe- 
sis. Hum. Brain Mapp. 26, 15-29. doi: 10.1002/hbm.20I13 

Greicius, M. D., Krasnow, B., Reiss, A. L., and Menon, V. (2003). Functional 
connectivity in the resting brain: a network analysis of the default mode 
hypothesis. Proc. Natl Acad. Sci. U.S.A. 100, 253-258. doi: 10.1073/pnas.01350 
58100 

Griswold, M. A., Jakob, P. M., Heidemann, R. M., Nittka, M., Jellus, V, 
Wang, J., et al. (2002). Generalized autocalibrating partially parallel acqui- 
sitions (GRAPPA). Magn. Reson. Med. 47, 1202-1210. doi: 10.1002/mrm. 
10171 

Hiifner, K., Stephan, T., Flanagin, V. L., Deutschlander, A., Stein, A., Kalla, 
R., et al. (2009). Differential effects of eyes open or closed in darkness on 
brain activation patterns in blind subjects. Neurosci. Lett. 466, 30-34. doi: 
10.1016/j.neulet.2009.09.010 

Hwang, K., Hallquist, M. N., and Luna, B. (2013). The development of hub archi- 
tecture in the human functional brain network. Cereb. Cortex 23, 2380-2393. 
doi: 10.1093/cercor/bhs227 

Jao, T, Vertes, R E., Alexander-Bloch, A. F, Tang, I. N., Yu, Y.-C, Chen, 
J.-H., et al. (2013). Volitional eyes opening perturbs brain dynamics and 
functional connectivity regardless of light input. Neuroimage 69, 21-34. doi: 
10.1016/j.neuroimage.2012. 12.007 



Frontiers in Human Neuroscience 



www.frontiersin.org 



August 2014 | Volume 8 | Article 645 | 7 



Liang et al. 



Sensorimotor fALFF of eyes conditions 



Marx, E., DeutscMander, A., Stephan, T., Dieterich, M., Wiesmann, M., 
and Brandt, T. (2004). Eyes open and eyes closed as rest conditions: 
impact on brain activation patterns. Neuroimage 21, 1818-1824. doi: 
10.1016/j.neuroimage.2003.12.026 

Marx, E., Stephan, T., Nolte, A., Deutschlander, A., Seelos, K. C, Dieterich, M., 
et al. (2003). Eye closure in darlmess animates sensory systems. Neuroimage 19, 
924-934. doi: 10.1016/S1053811903001502 

Mazard, A., Laou, L., Joliot, M., and Mellet, E. (2005). Neural impact of the seman- 
tic content of visual mental images and visual percepts. Brain Res. Cogn. Brain 
Res. 24, 423-435. doi: 10.1016/j.cogbrainres.2005.02.018 

Mazard, A., Mazoyer, B., Etard, O., Tzourio-Mazoyer, N., Kosslyn, S. M., 
and Mellet, E. (2002). Impact of fMRI acoustic noise on the functional 
anatomy of visual mental imagery. /. Cogn. Neurosci. 14, 172-186. doi: 
10.1162/089892902317236821 

Mcavoy, M., Larson-Prior, L., Ludwikow, M., Zhang, D., Snyder, A. Z., Gusnard, D. 
L., et al. (2012). Dissociated mean and functional connectivity BOLD signals in 
visual cortex during eyes closed and fixation. /. Neurophysiol. 108, 2363-2372. 
doi: 10.1152/jn.00900.2011 

Mcavoy, M., Larson-Prior, L., Nolan, T. S., Vaishnavi, S. N., Raichle, M. E., and 
D'Avossa, G. (2008). Resting States Affect Spontaneous BOLD Oscillations 
in Sensory and Paralimbic Cortex. /. Neurophysiol. 100, 922-931. doi: 
10.1152/jn.90426.2008 

Naselaris, T., Kay, K. N., Nishimoto, S., and Gallant, J. L. (2011). 
Encoding and decoding in fMRI. Neuroimage 56, 400—410. doi: 
I0.I016/j.neuroimage.20I0.07.073 

Norman, K. A., Polyn, S. M., Detre, G. 1., and Haxby, J. V. (2006). Beyond 
mind-reading: multi-voxel pattern analysis of fMRI data. Trends Cogn. Sci. 10, 
424-430. doi: 10.1016/j.tics.2006.07.005 

Poldrack, R. A., Halchenko, Y. O., and Hanson, S. 1. (2009). Decoding the large- 
scale structure of brain function by classifying mental States across individuals. 
Psychol. Sci. 20, 1364-1372. doi: 10.1111/j.l467-9280.2009.02460.x 

Qin, R M., Grimm, S., Duncan, N. W., Holland, G., Guo, ). S., Fan, Y., et al. (2013). 
Self-specific stimuli interact differently than non-self-specific stimuli with eyes- 
open versus eyes-closed spontaneous activity in auditory cortex. Front. Hum. 
Neurosci. 7:437. doi: 10.3389/fnhum.2013.00437 

Raichle, M. E., Macleod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and 
Shulman, G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. 
U.S.A. 98, 676-682. doi: 10.1073/pnas.98.2.676 

Shen, H., Wang, L., Liu, Y., and Hu, D. (2010). Discriminative analy- 
sis of resting-state functional coimectivity patterns of schizophrenia using 
low dimensional embedding of fMRI. Neuroimage 49, 3110-3121. doi: 
10.1016/j.neuroimage.2009.1 1.01 1 

Simmons, W. K., Avery, J. A., Barcalow, J. C., Bodurka, J., Drevets, W. C., and 
Bellgowan, P. (2013). Keeping the body in mind: Insula functional organization 
and functional connectivity integrate interoceptive, exteroceptive, and emo- 
tional awareness. Hum. Brain Mapp. 34, 2944-2958. doi: 10.1002/hbm.22113 

Song, X. W., Dong, Z. Y, Long, X. Y, Li, S. E, Zuo, X. N., Zhu, C. Z., et al. (2011). 
REST: a toolkit for resting-state functional magnetic resonance imaging data 
processing. PLoSOJV£6:e25031 . doi: 10.1371/journal.pone.0025031 

Uludag, K., Dubowitz, D. ]., Yoder, E. 1., Restom, K., Liu, T. T., and Buxton, R. 
B. (2004). Coupling of cerebral blood flow and oxygen consumption during 
physiological activation and deactivation measured with fMRI. Neuroimage 23, 
148-155. doi: 10.1016/j.neuroimage.2004.05.013 

Vergun, S., Deshpande, A. S., Meier, T. B., Song, )., Tudorascu, D. L., Nair, V. A., 
et al. (2013). Characterizing functional connectivity differences in aging adults 
using machine learning on resting state fMRI data. Front. Comput. Neurosci. 
7:38. doi: 10.3389/fncom.2013.00038 



Wang, J. H., Zuo, X. N., Gohel, S., Milham, M. R, Biswal, B. B., and He, Y. (2011). 
Graph theoretical analysis of functional brain networks: test-retest evaluation 
on short- and long-term resting-state functional MRI data. PLoS ONE 6:e21976. 
doi: 10.1371/iournal.pone.0021976 

Wang, K., Jiang, T., Yu, C, Tian, L., Li, J., Liu, Y., et al. (2008). Spontaneous activity 
associated with primary visual cortex: a resting-state FMRI study. Cereb. Cortex 
18, 697-704. doi: 10.1093/cercor/bhml05 

Wu, L., Eichele, T., and Calhoun, V. D. (2010). Reactivity of hemodynamic 
responses and functional connectivity to different states of alpha syn- 
chrony: a concurrent EEG-fMRI study. Neuroimage 52, 1252—1260. doi: 
10.101 6/j.neuroimage.20 1 0.05.053 

Xu, R, Huang, R., Wang, )., Van Dam, N. T., Xie, T., Dong, Z., et al. 
(2014). Different topological organization of human brain functional net- 
works with eyes open versus eyes closed. Neuroimage 90, 246-255. doi: 
10.1016/j.neuroimage.2013.12.060 

Yan, C, Liu, D., He, Y., Zou, Q., Zhu, C, Zuo, X., et al. (2009). Spontaneous 
brain activity in the default mode network is sensitive to different resting-state 
conditions with Umited cognitive load. PLoS ONE 4:e5743. doi: 10.1371/jour- 
nal.pone.0005743 

Yang, H., Long, X. Y., Yang, Y., Yan, H., Zhu, C. Z., Zhou, X. R, 
et al. (2007). Amplitude of low frequency fluctuation within visual areas 
revealed by resting-state functional MRI. Neuroimage 36, 144-152. doi: 
10. 10 1 6/j.neuroimage.2007.01 .054 

Zang, Y. P., He, Y., Zhu, C. Z., Cao, Q. J., Sui, M. Q., Liang, M., et al. (2007). 
Altered baseline brain activity in children with ADHD revealed by resting-state 
functional MRI. Brain Dev. 29, 83-91. doi: 10.1016/j.braindev.2006.07.002 

Zeng, L. L., Shen, H., Liu, L., Wang, L. B., Li, B. J., Fang, R, et al. (2012). Identifying 
major depression using whole-brain functional connectivity: a multivariate 
pattern analysis. Brain 135, 1498-1507. doi: 10.1093/brain/aws059 

Zhang, D., Liu, B., Chen, )., Peng, X., Liu, X., Fan, Y, et al. (2013). Determination of 
vascular dementia brain in distinct frequency bands with whole brain functional 
connectivity patterns. PLoS ONE 8:e54512. doi: 10.1371/journaLpone.0054512 

Zou, Q. H., Long, X. Y, Zuo, X. N., Yan, C. G., Zhu, C. Z., Yang, Y H., et al. 
(2009). Functional connectivity between the thalamus and visual cortex under 
eyes closed and eyes open conditions: a resting-state fMRI study. Hum. Brain 
Mapp 30, 3066-3078. doi: 10.1002/Hbm.20728 

Zou, Q. H., Zhu, C. Z., Yang, Y., Zuo, X. N., Long, X. Y., Cao, Q. J., el al. (2008). 
An improved approach to detection of amplitude of low-frequency fluctuation 
(ALFF) for resting-state fMRI: fi-actional ALFR /. Neurosci. Methods 172, 
137-141. doi: 10.1016/i.ineumeth.2008.04.012 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 14 March 2014; accepted: 02 August 2014; published online: 20 August 2014. 
Citation: Liang B, Zhang D, Wen X, Xu P, Peng X, Huang X, Liu M and Huang R 
(2014) Brain spontaneous fluctuations in sensorimotor regions were directly related to 
eyes open and eyes closed: evidences from a machine learning approach. Front. Hum. 
Neurosci. 8:645. doi: 10.3389/fnhum.2014.00645 

This article was submitted to the journal Frontiers in Human Neuroscience. 
Copyright © 2014 Liang, Zhang, Wen, Xu, Peng, Huang, Liu and Huang. This is an 
open-access article distributed under the terms of the Creative Commons Attribution 
License (CCBY). The use, distribution or reproduction in other forums is permitted, 
provided the original author(s) or licensor are credited and that the original publica- 
tion in this journal is cited, in accordance with accepted academic practice. No use, 
distribution or reproduction is permitted which does not comply with these terms. 



Frontiers in Human Neuroscience 



www.frontiersin.org 



August 2014 | Volume 8 | Article 645 | 8 



